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The discovery of the charge-density-wave formation in the high-Tc cuprate superconductors has activated in¬ 
tensive theoretical studies for the pseudogap states. However, the microscopic origin of the charge-density-wave 
state has been unknown so far since the many-body effects beyond the mean-field-level approximations, called 
the vertex corrections, are essential. Toward solving this problem, we employ the recently developed functional 
renormalization-group method, by which we can calculate the higher-order vertex corrections in a systematic 
and unbiased way with high numerical accuracy. We discover the critical development of the p-orbital-density- 
wave (p-ODW) instability in the strong-spin-fiuctuation region. The obtained p-ODW state possesses the key 
characteristics of the charge ordering pattern in Bi- and Y-based superconductors, such as the wave vector par¬ 
allel to the nearest Cu-Cu direction, and the Tsymmetry form factor with the antiphase correlation between px 
and py orbitals in the same unit cell. In addition, from the observation of the beautiful scaling relation between 
the spin susceptibility and the p-ODW susceptibility, we conclude that the main driving force of the density 
wave is the Aslamazov-Larkin vertex correction that becomes very singular near the magnetic quantum-critical 
point. 

PACS numbers: 74.72.Kf, 74.20.-z, 74.40.Kb, 75.25.Dk 


I. INTRODUCTION 

Understanding of the exotic electronic states in underdoped 
cuprate high-Tc superconductors has been one of the great¬ 
est challenges in condensed-matter physics. Especially, the 
charge-density-wave state observed in various compounds is 
recognized as a key feature of underdoped cuprate supercon¬ 
ductors [1-21]. The x-ray measurements have succeeded in 
the direct observation of the charge-density waves below ^ 
200 K [8, 9]. By means of the resonant x-ray scattering 
and STM measurements, it has now been recognized in many 
compounds that (i) the ordering vector of the charge-density 
wave is of axial type with the wave vector Qa = (^a^O) 
[8, 9, 11-15, 17-21], and (ii) the charge modulation emerges 
mainly on the oxygen p orbital [12, 17, 20, 21]. In addition, 
(iii) the symmetry of the charge-order pattern is ti-wave type, 
in which the modulations are antiphase between Px and Py 
orbitals in the intra unit cell [12, 17, 20, 21]. The charge- 
modulation pattern [12, 17, 20, 21] is depicted in Fig. 1(a), 
which can be regarded as the (i-symmetry p-orbital charge- 
density-wave (p-ODW) state. The density-wave wave vector 
Qa corresponds to wave vector connecting the neighboring 
“hot spots” shown in Fig. 1(b). The microscopic derivation 
of the charge-density-wave state with the key properties (i)- 
(iii) based on the realistic Hubbard model has been desired 
for years. 

Theoretical studies for the pseudogap phenomena have 
been performed by considering various strong-correlation ef¬ 
fects [22-26]. Motivated by the discovery of the density 
wave below ^ 200 K [8-21], many scenarios of the spin- 
fluctuation-driven nematic order have been proposed in Refs. 
[27-44], based on various single-orbital models. Especially, 
the bond-density-wave (BDW) state [28-37], the pair-density- 
wave (PDW) state [38-42], and the composite charge orders 
[43, 44] have been studied. The BDW state is deflned as the 
density wave with (i-symmetry form factor. As the driving 
force of the BDW state, the authors in Refs. [28-30, 33, 36] 


had focused on the Maki-Thompson vertex correction (VC), 
in analogy to the (i-wave superconductivity driven by the 
Maki-Thompson VC in the Eliashberg theory. The Maki- 
Thompson process for the BDW instability is shown in Fig. 
1(c), which is the flrst-order term with respect to the spin 
susceptibility. In general, the Maki-Thompson VC gives the 
density-wave instabilities at wave vectors q = Qa or Qd- 
The wave vectors Qa and Qd connect the hot spots on the 
Fermi surface (FS) shown in Fig. 1(b). The enhancement of 
the BDW susceptibility was supported by the renormalization- 
group (RG) method in the weak-coupling region [31, 37]. 
However, the predominant wave vector is of the diagonal 
type, Qd, inconsistently with experimentally-observed axial 
nematic order. It was also pointed out that the BDW instabil¬ 
ity driven by the Maki-Thompson process does not dominate 
over the superconducting instability [43, 45]. Also, the PDW 
state has been considered as the origin of the pseudogap phase 
[39, 40]. The PDW is formed by the linear combination of the 
Cooper pairs with flnite momenta. The PDW is also induced 
by the Maki-Thompson process [41,42], although the induced 
charge modulation has momentum 2Qa [40]. To realize the 
density wave with Qg,, several types of the composite charge 
order parameters have been proposed [41-44]. However, the 
d symmetry of its form factor has not been explained. 

Quite recently, the Aslamazov-Larkin VC was suggested to 
be more important near the magnetic quantum-critical point 
[46]. The Aslamazov-Larkin VC is the second-order term 
with respect to the spin susceptibility [Fig. 1(d)]. The pre¬ 
dominant wave vector driven by this VC is of the axial type, 
Qg,, due to the important scattering process shown in Fig. 1(b). 
However, only the single Aslamazov-Larkin process had been 
studied in Ref. [46]. In addition, the obtained form factor 
is given by a complex mixture between the d- and p-orbital 
charge densities. 

In order to settle down the controversy on the driving force 
of the density wave, we have to employ a sophisticated the- 
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FIG. 1. (a) Schematic charge distribution in the cZ-symmetry p-ODW 
state with the wave vector Qa = (O.Stt, 0). The size of the orbitals 
represents the charge density, (b) The Fermi surface of the cuprate 
superconductors. The density-wave instabilities due to the VCs 
emerge at Qd = (^d, (^d) and Qa = ((^a, 0), both of whieh connect 
the hot spots, (e) The Maki-Thompson (MT) process [30, 33, 36] 
and (d) the Aslamazov-Larkin (AL) process [46], which are the first- 
order and second-order terms with respect to the spin susceptibility, 
respeetively. The scattering process in the Maki-Thompson VC, ex¬ 
pressed by G{p) G{p + Qd) G{p') G{p' + Qd), is shown by green 
lines in (b). Also, the important scattering process in the Aslamazov- 
Larkin VC, which is expressed by G{k) G{k + Qa) G{k — q') in 
the three-point vertex, is shown by blue lines where q' ^ Qs (the 
wave vector of the spin fluctuation). 


oretical method to calculate higher-order diagrams, including 
both Maki-Thompson and Aslamazov-Larkin VCs, in a sys¬ 
tematic and unbiased way. For this purpose, the functional RG 
method would be the most appropriate theoretical technique. 
This method enables us to calculate various types of VCs up 
to the parquet-approximation level, in which the infinite series 
of the Maki-Thompson and Aslamazov-Larkin VCs are in¬ 
cluded. The RG method has been successfully applied in the 
one-dimensional electron systems [47, 48] and has been de¬ 
veloped as a powerful method for two-dimensional strongly- 
correlated electron systems [49-57]. 


In this paper, we employ the recently developed im¬ 
proved functional RG method, called the RG+constrained 
RPA (cRPA) method [58, 59], in order to tackle the unsolved 
theoretical problem on the charge-density-wave state. By uti¬ 
lizing this method, we can evaluate the momentum depen¬ 
dence of susceptibilities with high accuracy. We examine 
various charge and spin susceptibilities in the realistic three- 
orbital d-p Hubbard model [60, 61]. We discover that the sus¬ 
ceptibility of the (i-symmetry p-ODW state is critically en¬ 
hanced at the wave vector Qa in the strong-spin-fluctuation 
region. The obtained p-ODW state explains satisfactorily the 
experimental key features (i)-(iii) listed above. Thus, we pre¬ 
dict that the p-ODW with d symmetry [Fig. 1(a)] is the origin 
of the density wave in cuprate superconductors. The beautiful 
scaling between the spin and p-ODW susceptibilities means 
that the main driving force of the charge-density instability at 
Qa is the Aslamazov-Larkin VC, which is more singular than 
the Maki-Thompson VC near the magnetic quantum-critical 
point. Therefore, the p-ODW in cuprates originates from the 
strong interference between the spin and orbital fluctuations, 
which is described microscopically as the vertex corrections. 

This paper is organized as follows. In Sec. II, we intro¬ 
duce the three-orbital d-p model and explain the key idea of 
the RG+cRPA theory. In Sec. Ill, we show the numerical re¬ 
sults on the spin, p-ODW, and (i-orbital BDW susceptibilities. 
In the weak-spin-fluctuation region, the diagonal density sus¬ 
ceptibilities are moderately enhanced by the Maki-Thompson 
VC, consistently with the previous results [31, 33, 36, 37]. 
In the strong-spin-fluctuation region, in contrast, the axial p- 
ODW susceptibility is critically enhanced by the Aslamazov- 
Larkin VC, dominating over the BDW susceptibility. In Sec. 
IV, we show the beautiful scaling relation between the spin 
and charge susceptibilities, which means that the instability at 
q = Qa is driven by the Aslamazov-Larkin VC. We briefly 
discuss the effect of the Coulomb interaction for p orbital. 
Section V is devoted to conclusions. Details of the technical 
calculations are given in Appendixes, where the comparison 
between the conventional patch-RG and the RG+cRPA theo¬ 
ries are also made for the present three-orbital d-p model. 

II. d-p HUBBARD MODEL AND RG+CRPA METHOD 

We investigate a standard three-orbital d-p Hubbard model 
[60, 61], shown in Fig. 1(a), which has been analyzed for un¬ 
derstanding the charge-density-wave state on the oxygen p or¬ 
bital [46, 62, 63]. Its Hamiltonian is given by 

Hdp = y] cL Ho{k) + Vd y] nd,j,tnd,j 4 , (1) 

j 

where creation operator 

for the electron on 6 ^ 3 . 2 _^ 2 , p^., andp^ orbitals with wave vec¬ 
tor k and spin a. The (i-electron density operator is rid^j^o- = 
(it ^(ij,cr- We only consider the on-site Coulomb interaction 
for the d orbital Ua and the effect of intersite Coulomb inter¬ 
action is briefly discussed later. For the kinetic term Ho{k), 
we use the first-principles hopping integrals for La 2 Cu 04 in 
Ref. [64] (See the note in Ref. [65]). In addition, we introduce 
the third-nearest d-d hopping —0.1 eV to make the FS closer 
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FIG. 2. (a) The Fermi surface and (b) the band structure in the present three-orbital d-p model. The low-energy region (\E\ < 0.5 eV) are 
denoted by the shaded areas. The iVp-patch discretization for Np = 64 is shown in (a), whereas we set Np = 128 in the numerical study 
of the present paper, (c) The obtained spin susceptibility for the d orbital for Ud = 4.40 eV, T = 20 meV and Ao = 0.5 eV. The 

peak positions are q — Qs — and (tt, tt — (5s) where 6s ~ O.SItt. (d) The Ud dependencies of Xma5[= given by the 

RG-i-cRPA method and by the RPA. The initial spin susceptibility given by the cRPA is also shown. No “constrained VC” discussed in Ref. 
[58] is included in the present RG-i-cRPA study. 


to Y- and Bi-based cuprates by following Ref. [46]. The band 
filling is set to n = + np = 4.9, corresponding to the hole 

filling X = 0.1. The FS and the band structure are shown in 
Figs. 2(a) and 2(b), respectively. 

In underdoped cuprates, the spin, charge, and orbital de¬ 
grees of freedoms are strongly coupled by the electron cor¬ 
relation. This fact had prevented the explanation of the 
experimentally-observed nematic charge order so far. In 
order to analyze the competing strong fiuctuations in the 
two-dimensional systems accurately, we apply the RG+cRPA 
method developed in Refs. [58, 59], in which the q dependen¬ 
cies of susceptibilities can be obtained with high accuracy in 
comparison with the conventional patch-RG method. 

In the RG+cRPA method proposed in Ref. [58], we set the 
initial cutoff Aq smaller than the bandwidth [Fig. 2(b)]. The 
scattering processes of electrons having energies \E\ < Aq 
are integrated within the one-loop RG scheme based on the 
conventional A^-patch RG method [50-52, 56]. The division 
of lower-energy Brillouin zone (\E\ < Aq) into Np patches 
is illustrated in Fig. 2(a). The scattering processes involving 
higher-energy states with jE"! > Aq are evaluated by the con¬ 
strained RPA (cRPA) accurately using fine k meshes and in¬ 
corporated into the initial values. This treatment is based on 
the natural assumption that the vertex corrections become sig¬ 
nificant only in the low-energy regions |E| ^ Aq. Due to 
this treatment, low numerical accuracy for the higher-energy 


processes inherent in the conventional patch-RG method is 
greatly improved in the RG+cRPA method. Another advan¬ 
tage of the RG+cRPA theory for multi-orbital systems is that 
the inter-band processes (so-called van Vleck contributions) 
can be included. Therefore, the susceptibilities at low tem¬ 
peratures are obtained with high accuracy in the RG+cRPA 
method. 

In the present analysis we set Aq = 0.5 eV and Np = 
128. We have checked the case for Aq = 1.0 eV and find 
that the results are essentially independent of Aq. In the 
present RG+cRPA method, the numerical accuracy is suffi¬ 
ciently improved in comparison with the conventional patch- 
RG method. However, the main results of this paper are robust 
and reproduced qualitatively even in the conventional patch- 
RG method. The direct comparison between the numerical 
results of the RG+cRPA and those of the conventional patch- 
RG method are made in Appendix A. 

III. NUMERICAL RESULTS 
A. Spin susceptibility 

First, we analyze the spin susceptibility using the 
RG+cRPA method, by following the procedure explained in 
Refs. [58, 59]. Because of the (i-orbital Coulomb interaction 
Ud, the spin fiuctuations develop only on the d orbital. The 
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(i-orbital spin susceptibility per spin is given as 

VP“(g) = 1dT {S^{q, t) Sdi-q, 0)), (2) 

where Sd{q, r) is the spin operator for the d orbital. The mo¬ 
mentum dependence of the obtained is shown in Fig. 

2(c). The strong spin fluctuations are realized at the incom¬ 
mensurate wave vectors Qs = (tt — Ss^tt) and 
The obtained q dependence of is consistent with the 

neutron measurements. As increasing Ud, the maximum of 
the spin susceptibility, develops monotonically and di¬ 
verges at Ud = 4.5 eV), as shown in Fig. 2(d). The 

value of will increase to the first-principles value Ud ^ S 
eV by including the spin fluctuation-induced self-energy [ 66 ]. 

As seen from Figs. 2(c) and 2(d), the contributions from 
the cRPA are small for Aq = 0.5 eV but quite important 
(especially for the four-point vertex [58, 59]) in order to de¬ 
rive reliable results. In order to verify the validity of the 
RG-i-cRPA theory, we analyzed the same d-p model by us¬ 
ing the conventional patch-RG method established in litera¬ 
ture [50-52, 56] in Appendix A: It is confirmed that the es¬ 
sential results given by the RG-fcRPA method presented in 
the main text are qualitatively reproduced by the conventional 
patch-RG method, whereas the numerical accuracy is well im¬ 
proved in the RG-fcRPA method. 

B. p-ODW susceptibility and d-orbital BDW susceptibility 

In contrast to spin susceptibility, any charge susceptibilities 
are not enhanced by Ud in the mean-field-level approxima¬ 
tions. Nonetheless, we reveal that the strong p-ODW insta¬ 
bility emerges in the present RG analysis, thanks to the VCs 
that are dropped in the RPA. Since Up = 0, the effective in¬ 
teraction on p orbitals that causes the p-ODW instability is 
originated from the many-body effects on the d orbital. The 
p-ODW susceptibility per spin is defined as 

Xa/3 (?) = 2 dr {na{q,T)ni3{-q,0)), (3) 

where (3 = x,y represent the Px^Py orbitals. According to 
the experimental analyses [12, 17], we introduce the suscepti¬ 
bilities for the p-ODW with d and s' symmetries: 

x7"\q) ^ xr”"(g)+xrr"(9) - 2xrr"(9), (4) 
x:;”"(g)=xr”"(9)+ xi7\q)+ 2xrr"(g)- (5) 

The susceptibility measures the development of the 

antiphase correlation between px and Py orbitals in the same 
unit cell. Aside from the conventional charge/orbital orders, it 
was pointed out in literature [28-30, 33, 35, 63] that the bond- 
density order with d-wave form factor is expected to develop 
in the Hubbard model. Therefore, we calculate the d-orbital 
BDW susceptibility: 

^BDW(q) ^ dr {B{q,T) B{-q,0)), (6) 

Biq) = X] /(fc + g/2) (7) 

k,a 


where f{k) is the d-wave form factor: f{k) = cos{kx) — 
cos{ky). The order parameter {B{q)) / 0 represents the 
bond-ordered state, which is equivalent to the modulation of 
the hopping integrals [35]. In the case of q = 0, the suscep¬ 
tibility X^^^ {q = 0) measures the nematic or Pomeranchuk 
instability [50-52, 56]. We also evaluate the conventional d- 
orbital charge susceptibility 

1 

^d-orb(g) = - dr {nd{q, r) nd(-q, 0 )), ( 8 ) 

where nd{q, r) is the density operator for the d orbital. Al¬ 
though the upper limit of ^d-orb 

(q) is 1/Ud in the present 
Hubbard model, x^~^^^{q) the one-loop RG may diverge 
unphysically if we set Ud too large. Thus, we always ver¬ 
ify the non-singular behavior of x^~^^^{q) to ensure that the 
adopted parameter value is within a valid range. We applied 
the finite-temperature RG formalism based on the sharp band¬ 
width cutoff and solved the RG equations down to A/ = 0. 
In the actual calculation of the RG equations, we introduce a 
lower-energy cutoff A/ = ttT for the four-point vertex func¬ 
tion in order to obtain stable numerical results [31]. The phys¬ 
ical meaning of this lower-energy cutoff is, for example, the 
suppression of the Cooper channel due to the impurity scatter¬ 
ing or the magnetic field. Experimentally, the charge-density- 
wave state is strongly stabilized by applying the magnetic field 
beyond 15 T [67, 68 ]. This fact means that the Cooper insta¬ 
bility is less important for the density-wave-formation mech¬ 
anism. 

In Fig. 3, we present the obtained p-ODW susceptibilities 
with d and s' symmetries, together with the d-orbital BDW 
susceptibility for Ud = 4.40 eV in (a)-(d), and Ud = 4.50 eV 
in (e)-(h). In the case of Ud = 4.40 eV, in which the spin sus¬ 
ceptibility is moderate (XmS ~ eV“^), both xS"°^^ 

xBDW are enlarged compared to the RPA results as shown 
in Figs. 3(a) and 3(c), whereas Xs"°^^ Fig- 3(b) is not en¬ 
hanced at all. Therefore, the p-ODW and BDW susceptibil¬ 
ities are moderately enhanced by the VCs that are neglected 
in the RPA. However, the highest peaks of both xS"°^^ 
xBDW are located at q = Qd, inconsistently with the axial 
nematic density wave in cuprates. 

In the case of Ud = 4.50 eV, in which the spin suscep¬ 
tibility is large (Xmix ~ 200 eV“^), both the p-ODW and 
d-orbital BDW susceptibilities shown in Figs. 3(e) and 3(g) 
possess large sharp peaks at q = Qa and Qd, which origi¬ 
nate from the VCs generated in the renormalization procedure. 
The most divergent density-wave susceptibility is 
q = Qa- In contrast, xl~^^^{q) is seldom enhanced. There¬ 
fore, the axial p-ODW with d symmetry shown in Fig. 1(a) 
is realized. We have verified that |Qa| = increases with 
hole doping in the present RG study, consistently with experi¬ 
ments in the Y-, Bi-, and Hg-based compounds. According to 
neutron inelastic scattering studies, Xmi5 is as large as ^ 500 
eV“^ in the slightly underdoped YBCO [69]. Therefore, the 
nematic density wave in cuprates is realized in the strong-spin 
fiuctuation region experimentally. 

We note that both xS"°^^ xB^^ have sub-dominant 
broad peaks at q = 0, indicating that “the q = 0 Pomer- 
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FIG. 3. The p-ODW and and d-orbital BDW [x^^^(^)] susceptibilities in the three-orbital Hubbard model with 

Ud = 4.40 eV in (a)-(d), and = 4.50 eV in (e)-(h). The temperature is T = 20 meV. The maximum of the spin susceptibility is XmS ~ 12 
eV“^ (200 eV“^) for Ud = 4.40 eV (4.50 eV). The RPA results are also shown for comparison. The axial wave vector is Qa ~ (O.STtt, 0) 
and the diagonal wave vector is Qd ~ (0.4l7r, 0.4l7r). Both peak positions Qa and Qd correspond to the wave vector connecting the hot 
spots shown in Fig. 1(b). The wave vectors Qa and Qd obtained from the RG-i-cRPA method change with carrier doping, by following the 
change in the vectors connecting the hot spots. 


anchuk instabilities” are also enhanced in the present model. 
This instability had been reported in the previous theoreti¬ 
cal studies [31, 70-72], and also observed experimentally in 
cuprates as the enhancement of the Big channel Raman re¬ 
sponse [70-72]. The temperature-flow RG scheme [53] would 
also be useful for the study of the Pomeranchuk instability 
[50, 54]. 


IV. DISCUSSIONS 

A. Scaling relation between p-ODW susceptibility 
and spin susceptibility 

In the inset of Fig. 4(a), we show the Ud dependencies of 
p-ODW and (i-orbital BDW susceptibilities at q = Qa- The 
p-ODW susceptibility exceeds the BDW one with increasing 
Ud. In the main flgure of Fig. 4(a), we plot both the p-ODW 
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(a) 



(b) 



FIG. 4. (a) The p-ODW and (i-orbital BDW susceptibilities, and at peak positions as functions of The short- 

dashed lines represent the RPA results, which are almost independent of Ud- starts to increase when XmS exceeds 10 eV“^ 

(shown by the arrow). The corresponding Stoner enhancement is XmS/xo ~ 40. The gradient of the dotted line is 1, so a beautiful scaling 
relation ^ XmaS is satisfied in the strong-spin-fiuctuation region (XmS > 10)- The inset shows the Ud dependencies of 

and at qr = (b) The d-symmetry and s'-symmetry p-ODW susceptibilities, Xd'^’^ {o) and {o)^ as well as the (i-orbital 

BDW and charge susceptibilities, and at qr = Qa, as functions of XmS- In the strong-spin-fiuctuation region, 

exhibits a critical enhancement, whereas the (i-orbital charge susceptibility remains small. Thus, the enhancement of Xd'° (<3a) is 

obtained within the reliable parameter range of the RG method. 


and (i-orbital BDW susceptibilities at gr = Qa and Qd as func¬ 
tions of Xmi5’ order to reveal the correlation between the 
spin and density susceptibilities. In the weak-spin-fluctuation 
region (Xmi5 ^ lb eV“^), both the p-ODW and (i-orbital 
BDW susceptibilities increase moderately. In the strong-spin- 
fluctuation region, only starts to increase drastically 

in proportion to XmiS- 1^ Ibis region, the highest peak of 
shifts to gr = consistently with the experimen¬ 
tal wave vector. The relation Xd~°^^(Qa) ^ xBDW(Qa) is 
robust against the choice of model parameters. 

The most important finding in Fig. 4(a) is that 
at gr = Qa well scales to Xmix Ib^ strong-fluctuation re¬ 
gion (Xmix ^ lb eV“^). This beautiful scaling relation is 
obtained in the wide range of model parameters. This fact 
indicates that the p-ODW is driven by the Aslamazov-Larkin 
VC that describes the strong interference between spin and 
orbital fluctuations. Also, as shown in Fig. 4(b), 
well as the conventional charge susceptibility for the d orbital 
X^-orb(gr) decrease with Xmix Ibe weak-spin-fluctuation re¬ 
gion, whereas they turn to increase slightly in the strong-spin- 
fluctuation region. The obtained relation < ^/Ud 

supports the reliability of the present RG result even for 
~ 200 eV-i. 

B. Why is the axial p-ODW susceptibility enlarged in the 
strong-spin-fluctuation region? 

In order to understand the physical origin of the p-ODW 
and, in addition, in order to confirm the validity of the present 
RG analysis, we also perform the diagrammatic analysis for 
the d-p Hubbard model. Hereafter, we demonstrate that 


the characteristic behaviors of the p-ODW and BDW sus¬ 
ceptibilities are well understood by taking VCs for the irre¬ 
ducible susceptibilities. The lowest-order Maki-Thompson 
and Aslamazov-Larkin VCs, Xmt{q) and Xal{q), are 
shown diagrammaticahy in Fig. 5(a). Their exact expressions 
are given in Appendix C. In the strong-spin-fluctuation re¬ 
gion, the Aslamazov-Larkin VC is scaled as [46] 

XAL(Qa) (X + OC OC xX, 

Q 

(9) 

where ^ is the magnetic correlation length, while the Maki- 
Thompson VC is scaled as 

XMT(<5a,d) OC + Qa.d) Oclll OC InX^^- 

( 10 ) 

Thus, the Aslamazov-Larkin VC is expected to dominate over 
the Maki-Thompson VC in the strong-spin-fluctuation region. 

Figure 5(b) shows the q dependencies of the Maki- 
Thompson and Aslamazov-Larkin VCs for the p-ODW sus¬ 
ceptibility, at Ud = 4.06 eV and T = 50 meV. In this case, 
the system is in the intermediate-spin-fluctuation region with 
Xmix ~ 24 eV“^. Here the Aslamazov-Larkin VC possesses 
the highest peak at gr = and the Maki-Thompson VC has 
the second highest peak at gr = Qd- Thus, the irreducible 
susceptibility 

(9) + (11) 

has the largest peak at gr = Qa due to the Aslamazov-Larkin 
VC. 
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FIG. 5. (a) The lowest-order Maki-Thompson VC XyiTiq) and Aslamazov-Larkin VC XAi.{q)^ Their exact mathematical expressions are 
given in Appendix C. (b) The momentum dependencies of the Maki-Thompson and Aslamazov-Larkin VCs for the p-ODW susceptibility, 
at T = 50 meV The magnitude of Ud is set to satisfy the Stoner factor as = 0.99. The noninteracting susceptibility 
shown, (c) The Maki-Thompson and Aslamazov-Larkin VCs for the cZ-orbital BDW susceptibility, for as = 0.99 with T = 50 meV. At 
q — Qd, the Maki-Thompson and Aslamazov-Larkin VCs are comparable. The inset shows that does not scale to XmaS- (d) 

The full susceptibilities and as functions of XmS- The scaling relation Xd'°’^^(Qa) oc XmS is well satisfied in the 

strong-spin-fluctuation region. For comparison, the RPA result for is also shown. 


We also show the Maki-Thompson and Aslamazov-Larkin 
VCs for the (i-orbital BDW susceptibility in Fig. 5(c). Since 
both VCs have large peaks at gr = Qd. the irreducible suscep¬ 
tibility 

+ XBDW(q) + (12) 


takes the highest peak at gr = Qd- As shown in the inset of 
Fig. 5(c), T>^^^(gr) at gr = Qd does not scale to Xmix’ which 
is consistent with the smallness of X^^"^ {q) i^ Fig- 4(a). 

Next, we derive the full susceptibilities from the irreducible 
susceptibilities. Figure 5(d) shows the obtained 
and at q = Qg, as functions of by changing 

Ud. Since Up = 0 in the present model, the p-ODW suscep¬ 
tibility is well approximated as ^ The 

obtained Xd~°^^(Qa) behaves very similarly to the RG result 
in Fig. 4(a). Especially, the scaling relation xr (<5a) oc 
Xmix is well reproduced by the diagrammatic analysis, due 
to large contribution from the Aslamazov-Larkin VC. In con¬ 
trast, the (i-orbital charge susceptibility 






(13) 


at gr = Qa shows a minimum at finite XmS’ similarly to 
the RG result in Fig. 4(b). This behavior is also understood 
analytically: In the weak-spin-fluctuation region, in which 


the VCs are negligible and therefore T>^‘°^^(gr) « Xo~^^^{q) 
is satisfied, the (i-orbital susceptibility decreases with Ud 
in proportion to 1/[1 + UdXo~^^^{q)]- 1^ Ihe strong-spin- 
fluctuation region, T>^“°^^(gr) increases drastically because of 
the Aslamazov-Larkin VC, and therefore x^~^^^{q) increases 
toward 1/f/d- 

Thus, we revealed that the characteristic behaviors of 
the p-ODW susceptibility in the present RG study, such 
as the peak position at gr = Qa and the scaling relation 
x2"°^^(Qa) Xmi5’ qualitatively understood by includ¬ 
ing the Aslamazov-Larkin VC into the RPA. This result is 
never trivial in that the higher-order VCs, unrestricted to the 
Maki-Thompson and Aslamazov-Larkin VCs, are systemati¬ 
cally produced in the RG theory. For example, the higher- 
order Maki-Thompson and Aslamazov-Larkin processes are 
included in the RG. Also, the spin and charge fluctuations and 
the four-point VCs are calculated consistently. Thus, the dom¬ 
inant role of the lowest-order Aslamazov-Larkin VC shown in 
Fig. 5 is confirmed in the present RG theory. 

Of course, the lowest-order Aslamazov-Larkin VC study 
cannot explain the RG results in many parts. For example, if 
only the lowest-order Aslamazov-Larkin VC is included, the 
relation X^^^(Qa) Xmix realized although it remains 
small in the RG results. Also, the (i-symmetry form factor 
in the p-ODW cannot simply be obtained by the lowest-order 
Aslamazov-Larkin VC. These facts indicate the importance of 
the higher-order diagrams included in the RG. 
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C. Effect of the Coulomb interactions for O sites on the 
p-ODW susceptibility 

The scaling relation Xd~^^^iQa) ^ XmiS obtained in Fig. 
4 indicates that the spin-density-wave order and the p-ODW 
order emerge simultaneously. However, the p-ODW will be 
realized in the paramagnetic state in the presence of small but 
finite Coulomb interaction between the nearest-neighbor p-p 
orbitals, Vpp. The p-ODW susceptibility is given as 


r'°'^{q-,Up,Vpp) 


l + Viq) xP-”b(q)’ 


(14) 


no. 6. The p-ODW susceptibility xV'^iQs^', Up, Vpp) [Eq. (14)] 
as a function of XmS for Cp = 0.7 eV. As an approximation, the 
irreducible susceptibility is given by the RG result for Up = V^p = 

0. 

Note that there is another type of the Aslamazov-Larkin VC 
given by the product of spin and charge (or orbital) propaga¬ 
tors ^ ^charge gives rise to the correc¬ 

tion to the spin susceptibility and is also included in the 
present RG formalism. From the data of spin susceptibility 
[Fig. 2(d)], it can be confirmed that such a VC does not have 
a strong effect. 


where represents the irreducible susceptibility with 

respect to Up (the p-orbital onsite interaction) and Vpp. The 
non-zero matrix elements of V{q) are Vxx = Vyy = Up and 
Vxy = Vyx = 8VppCos(ga:/2)cos(g^/2). When Up, Vpp are 
very small, x^"°^^(g^) is safely approximated by the RG sus¬ 
ceptibility for Up = Vpp = 0. Figure 6 shows the obtained 

xV'^iQv, Up, Vpp) for Vpp = 0.2 and 0.3 eV with Up = 0.7 
eV, which diverges even when finite. Now, we explain 

why the tiny Vpp critically enhances the p-ODW susceptibil¬ 
ity. Considering the relation Xx~x^^Xl~y^^^ ~ [Xx~y^^]‘^ 
tained in the present RG at gr = Qa. the d-symmetry p-ODW 
susceptibility is approximately given as 


x7"''iQv,Up,Vpp) 


xr\Q^) 


^ + Up [x'xUiQa) + X^UiQa)] 


2)4,(Qa)xSr’’(a 


(15) 


which is enhanced by Vpp since is negative. The 

enhancement due to Vpp dominates over the suppression due 
to Up, if Up ^ SVpp cos (5a. In real materials, finite electron- 
phonon coupling would also enlarge the p-ODW suscepti¬ 
bility in the strong-spin-fluctuation region. Therefore, the 
Aslamazov-Larkin VC accounts for the p-ODW ordering in 
the pseudogap region of the cuprate superconductors. 

V. CONCLUSIONS 

In this paper we applied the RG-fcRPA theory to the 
three-orbital d-p Hubbard model and discovered that the d- 
symmetry p-ODW susceptibility critically develops at q = 
Qa in the strong-spin-fiuctuation region. The main result 
is shown in Figs. 3(e)—3(h). The obtained p-ODW state 
has the characteristics: (i) the ordering vector is axial type 
Qa = (^a: 0). (ii) charge modulation occurs on the oxygen 
p orbital, and (iii) the symmetry of the p-orbital order pattern 
is of the (i-wave type. Therefore, the present RG-i-cRPA the¬ 
ory reproduced satisfactorily the experimental charge-density- 
wave state. The p-ODW originates from the strong interfer¬ 
ence between the spin and orbital susceptibilities. Such an 
interference is the main characteristics of the electronic states 
in underdoped cuprates. 


In the previous scenarios of the spin-fiuctuation-driven 
density-wave states, such as BDW [28-37], PDW [38-42], 
and the composite orders [43, 44], the Maki-Thompson VC 
has been studied as possible origins [28-30, 33, 36, 41-43]. 
In the weak-spin-fiuctuation region, the Maki-Thompson VC 
is dominant and the BDW instability with q = is ob¬ 
tained, consistently with Refs. [28-31, 33, 36, 37]. How¬ 
ever, in the strong-spin-fluctuation region, the p-ODW with 
the axial wave vector Qa is critically enhanced: The obtained 
beautiful scaling between the spin and p-ODW susceptibili¬ 
ties means the important role of the Aslamazov-Larkin VC. 
Therefore, the p-ODW in cuprates originates from the strong 
interference between the spin and orbital fiuctuations, which 
is microscopically described as the vertex corrections. 

Despite that the p-ODW and the (i-orbital BDW have the 
same d symmetry, we found that the p-ODW susceptibility is 
strongly enhanced while the enhancement of the BDW sus¬ 
ceptibility is moderate, as seen from Fig. 4(a). This fact im¬ 
plies that the charge-density modulation in cuprates is not due 
to the bond modulation of the d orbital, but due to the charge 
modulation on the p orbital. In fact, the recent x-ray diffrac¬ 
tion study reported that the sizable oxygen-site displacements 
occur in the charge-density-wave state, whereas the Cu-site 
displacements are very small [73]. This result supports the 
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p-ODW scenario proposed in this paper. 

The RG+cRPA theory developed in this study will be use¬ 
ful for analyzing unsolved problems in strongly-correlated 
electron systems. For example, it is interesting to apply the 
RG-fcRPA theory for Fe-based superconductors in order to 
understand the origin of the electronic nematic state [74-77]. 
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FIG. 8. The RG equations for the three-point vertex R{q; 0, 0'), 
where 0 is the patch index and I is the scaling parameter. The slashed 
(crossed) line represents an electron propagation having energy on 
Ai-^di < \E\ < Ai {\E\ > Ai) where Ai = Aoe~K The patch 
indices are determined so that the momenta hi and hi + q {i — 1^2) 
are respectively on the and 0- patches. Therefore, ki ^ p(0O 
and ki-\-q / p(0-). 


Appendix A: COMPARISON BETWEEN RG+cRPA 
AND THE CONVENTIONAL PATCH-RG SCHEME 

In the main text, we have analyzed the susceptibilities us¬ 
ing the recently developed RG-fcRPA theory [58, 59]. In the 
conventional patch-RG method, the numerical accuracy for 
the higher-energy processes becomes worse because of the 
large patch radius [50-52, 56]. In order to improve the nu¬ 
merical accuracy, in the RG-i-cRPA method, the higher-energy 
processes are calculated accurately within the cRPA by in¬ 
troducing fine k meshes. For this reason, the q dependen¬ 
cies of susceptibilities are obtained very accurately. Although 
the effect of VCs is underestimated in the RG-i-cRPA method, 
this underestimation is not serious since the VCs are impor¬ 
tant mainly in the lower-energy processes. In this section, we 
make a direct comparison between the numerical results of the 
RG-i-cRPA and those of the conventional patch-RG method in 
order to confirm the validity and reliability of the RG-fcRPA 
theory. 

The band structure and the patch discretization in the con¬ 
ventional patch-RG method are shown in Figs. 7(a) and 7(b), 
respectively. The Brillouin zone is divided into the patch seg¬ 
ments with respect to the angle 4>i [Fig. 7(b)]. Here, the radial 
variable from the center of each patch, p{(t)i), is neglected 
in the vertices, that is, r({/ci}) = T{{kr^i,(t)i}) 


(a) (b) 



FIG. 7. (a) Band structure and (b) Fermi surface in the present model. 
The discretization of the Brillouin zone {Np — 64) in the conven¬ 
tional patch-RG scheme are shown. In each patch, the center mo¬ 
mentum p{(j)i) is shown as the open circle. The initial cutoff Aq is 
set to the bandwidth of the conduction band, and the contributions 
from the valence bands are neglected. 


for the four-point vertex and R{q;ki,k 2 ) ^ , ^ 2 ) for 

the three-point vertex [50-52, 56]. The justification of this 
approximation is frequently ascribed to the simple scaling ar¬ 
gument in the weak-coupling limit. However, the k^ depen¬ 
dencies are quantitatively important in reality, since the “mo¬ 
mentum mismatch” problem is serious in the higher-energy 
processes Moreover, this problem becomes more serious in 
multi-orbital systems, since the orbital components of each 
band, called the “orbital makeup” [56], depend on /cr in usual 
multi-orbital models. 

The conventional patch-RG scheme can be reproduced in 
the present RG-bcRPA method by neglecting the valence-band 
contributions and by setting the cutoff Aq to the bandwidth 
of the conduction band [see Fig. 7(a)]. The RG equation for 
R{q] 01, ^ 2 ) is shown diagrammatically in Fig. 8. The sus- 
ceptibilities X®P“(q), and 

^d-orb 

{q) obtained by the conventional patch-RG method for 
Ud = 9.5 eV are shown in Fig. 9. As shown in Fig. 9(a), 
strong antiferromagnetic fluctuations are obtained, although 
the obtained incommensurate peak position of is at 

q = (tt — (5s, TT — (5s), differently from the experimental peak 
positions aX q = — 6s) and (tt — 6s, tt) that are repro¬ 

duced by the RPA as well as the RG-i-cRPA theory [see Fig. 
2(c)]. 

Figure 10 shows the inverse of the spin susceptibilities 
[xspin(g)]^G+cRPA and [x''P'"(q)]conv-RG, obtained in the 
RG-fcRPA and conventional patch-RG methods, respectively. 
In the case of Ud ^ 0, [x®^^^(q)]rg+cRPA almost perfectly 
reproduces the RPA susceptibility Xrpa( 9^) = x°(9)/(l - 
Ux^{q)), where x^{q) is the bare susceptibility. It is note¬ 
worthy that the result of the RG-i-cRPA method for Aq = 0.5 
eV is very similar to that for Aq = 1.0 eV. In contrast, 
[X^P'"(Q)]conv -RG gives under-estimated value since the van 
Vleck contribution is dropped. In addition, [x®^^^(9^)]conv-RG 
also deviates from the RPA restricted to the conduction- 
band contribution, Xrpa cond(9^)- origin of this de¬ 

viation is that the orbital components, called the “orbital 
makeup”[56, 78], largely depend on k^ because of the large 
patch radius in the conventional patch-RG method. Thus, the 
RG-i-cRPA method is superior to the conventional patch-RG 
method in the numerically accuracy for Ud ^ 0. 

Next, we compare the results of two RG methods for fi- 
nite Ud- As shown in Fig. 10, each l/[x®'’‘”(q)]RG+cRPA and 
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(0,0) (jt,0) (jt,jt) (0,0) 

no. 9. (a) (b) xV'iq), (c) (d) xr'^(9). 

and (e) obtained by the conventional patch-RG method. 

We SQi Ud = 9.5 eV, T = 20 meV, and Np = 64. For compari¬ 
son, the RPA result for (q) with = 4 eV is also shown in 
(a), and the RPA results for Ud = 9.5 eV are shown in (b)-(e). As 
shown in (b) and (c), c?-symmetry p-ODW instability at qr = Qa de¬ 
velops strongly. These results are consistent with the results of the 
RG-i-cRPA method explained in the main text. 


l/[VP“(9)]conv -RG is concave as functions of Ud because of 
the VC generated in the three- and four-point vertices. We 
see that Uf « 4.5 — 5 eV in the RG+cRPA method, whereas 
« 9.5 eV in the conventional patch-RG method. To un¬ 
derstand such large difference in in Fig. 8, we comment 
on the “momentum mismatch” between the pair of the Green 
functions G{k 2 )G{k 2 + q) and F, R, originating from the 
fact that the momenta at the centers of the patches 02 and 



FIG. 10. The Ud dependence of XmS obtained by the conventional 
patch-RG method, RG-i-cRPA methods, and RPA. In the RG-i-cRPA, 
the results with the different choices of parameters are shown. In 
the weak-coupling limit, the RG-i-cRPA results well reproduce the 
RPA results, irrespective to the choice of Np. In addition, the Ao 
dependence of the results is weak in the RG-i-cRPA. 


02 ’ P( 02 ) and p( 02 ), are different from k 2 and k 2 + q, re¬ 
spectively. This mismatch leads to violation of the momen¬ 
tum conservation at each bare vertex Ud in the Feynman di¬ 
agrams. This defect is not serious in the RG+cRPA method 
with Ao <C VFband- Considering that the self-energy cor¬ 
rection is dropped in the present RG methods, large U^^ in 
the conventional patch-RG method might partially originate 
from the momentum mismatch. On the other hand, U^^ in the 
RG+cRPA method is underestimated since the vertex correc¬ 
tions with higher-energy processes (\E\ > Aq) are included 
only partially. In spite of this underestimation, large p-ODW 
susceptibility is induced by the VC in the RG+cRPA method. 

We stress that the results of the RG+cRPA method and 
those of the conventional patch-RG method are qualitatively 
equivalent, nonetheless of the large difference in Uf. In the 
RG+cRPA method, the results are almost unchanged in the 
both cases of Aq = 0.5 eV and Aq = I.O eV In addition, 
the RG+cRPA results with Np = 64 and 128 are almost iden¬ 
tical, indicating that the the numerical convergence has been 
achieved at A/p « 64. We note that various RG methods have 
been proposed to resolve the above-mentioned problems in 
the conventional patch-RG method, such as the singular-mode 
functional renormalization group [78]. 

Regarding the charge susceptibilities, both and 

xBDW(^gf) possess large sharp peaks at q = Qa and 
as shown in Figs. 9(b) and 9(c) dXUd = 9.5 eV This be¬ 
havior is essentially identical to that obtained by RG+cRPA 
method (Fig. 3). This result supports the validity of the present 
RG+cRPA theory. Since the van Vleck contributions are 
dropped in the conventional patch-RG method, both 
and in Fig*- 9(b) and 9(d) are much smaller than the 

RPA results except near the peak positions. In the RG+cRPA, 
in contrast, the RG+cRPA results coincide well with the RPA 
results if q is away from the peak positions, as shown in Fig. 
3. 

Figure II shows the scaling relation between the p- 
ODW susceptibility and Xmi5 obtained from the conventional 
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FIG. 11. Peak values of and as functions of 

XmS obtained from the conventional patch-RG method. The gradi¬ 
ent of the dotted line is 1, so beautiful scaling relation ^ 

Xma5 is satisfied in the strong-coupling regime. The inset shows the 
Ud dependencies of and at qr = Qa- 


patch-RG method. The obtained relations oc 

Xmix as well as x^a^iQa) » in the strong- 

Spin-fluctuation region are qualitatively consistent with the 
RG+cRPA results shown in Fig. 5. Therefore, the develop¬ 
ment of the (i-symmetry axial p-ODW instability due to the 
spin fluctuation, which is the main message of the main text 
derived from the RG+cRPA, is confirmed also by the conven¬ 
tional patch-RG method. 

In summary, the results given by the RG+cRPA method 
presented in the main text are qualitatively reproduced by the 
conventional patch-RG method. Therefore, the validity of the 
RG+cRPA method is verified by the conventional patch-RG 
method, the validity of which has been confirmed in litera¬ 
ture [50-52, 56]. The numerical accuracy is well improved 
in the RG+cRPA method, as recognized by the coincidence 
with the RPA result in the weak-coupling region shown in Fig. 
10. This improvement is achieved by calculating the higher- 
energy processes accurately within the cRPA by introducing 
the fine k meshes. As seen from Figs. 2(c) and 2(d), the con¬ 
tributions from the cRPA are small for Aq = 0.5 eV but quite 
important (especially for the four-point vertex [58, 59]) in or¬ 
der to derive reliable results from the RG equations. 

Appendix B: SCALING FLOWS OF EFFECTIVE 
INTERACTIONS 

In this section, we examine the effective interactions for p- 
ODW and ti-orbital BDW. The effective interaction for the d- 
orbital BDW state is given by taking into account the (i-form 
factor 

rBDW(g) = (fe, fe'; k + q,k'-q) 

k,k' 

X u*^{k) ul{k') Ud{k + q) Ud{k' - q) 
x/(fe + q/2)/(fe'-q/2), (Bl) 



FIG. 12. The effective interactions (Bl) and (B2) at the scaling 
parameter / = 0, 2, and 3, with Ud = 4.50 eV, T = 20 meV, 
Np = 128. In this case, Xd~°’^^(Qa) is larger than x^^^(Qd) [see 
Fig. 3(h)]. 


where V^(k^k']k + q^k' — q) represents the renormalized 
charge vertex. This expression is essentially the same as ex¬ 
amined in the single-orbital case [31]. Here, U(i{k) represents 
the unitary transformation from the conduction band basis to 
the (i-orbital basis. The summations over k and k' are re¬ 
stricted to the low-energy scattering processes, in which all 
the momenta {k^k' ^k^ q^k' — q) are near the Fermi surface. 
This interaction becomes zero if the bare Hubbard interaction 
Ud is used for F^ [31]. In the present case, be¬ 

comes small but nonzero due to the cRPA contributions. 

In the similar manner, the effective interaction for the p- 
ODW state can be introduced as 

-KvxyiQ)-'^'yxyx(Q)^ (B2) 

k,k' 

X tx* (/c) Up{k') u^{k -b q) us{k' - q), (B3) 

where Ua{k) (a = x,y) represents the unitary transformation 
from the band basis to the Px- or -orbital basis. 

As has been discussed in Ref. [31], the negative values of 
the effective interactions indicate the precursor of the corre¬ 
sponding instability. Actually, if the effective interaction be¬ 
comes negative, the rhs of the RG equations for the three-point 
vertex [58] becomes positive and then the three-point vertex 
can increase. Once the three-point vertex is enhanced, the cor¬ 
responding susceptibility is also enhanced [58]. 

In Fig. 12, we show the development of the effective inter¬ 
actions for several scaling parameters I = ln(Ao/A) in the 
case ofU = 4.50 eV. At the initial point (I = 0), these ef¬ 
fective interactions given by the cRPA are small and almost 
independent of q. After the renormalization (/ = 2,3), the ef¬ 
fective interactions at Qa and Qd become negative, implying 
the development of the susceptibilities. At I = 3, the inter¬ 
action at q = Qa is stronger than that at q = Qd- We And 
that the overall profile of r^^^((3a) is similar to r^"°^^(Qa), 
except for the constant factor. 

In contrast, in the weak-fluctuation case (Ud = 4.40 eV), 
the interaction at q = Qa becomes weaker than that at q = 
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MT-VC 


{k+qliy 


AL-VC 



k'+qH) 


u^ip+q) 


uhip+q) 




FIG. 13. Detailed diagrammatic expressions of the Maki-Thompson and Aslamazov-Larkin VCs. Their analytic equations are shown in Eqs. 
(C3)-(C6). 


Qd- This relation of the effective interactions is in accordance 
with that of the peak positions of the susceptibility. 

Thus, we understand that the enhancements of 
and ^ = Qa, Qd, and 0 originate from the large 

negative and after the renormalization. 

However, the q dependencies of and 

not consistent with those of the susceptibilities. For example, 
the peak position of each and is Q = Qa 

or Qd for Ud = 4.50 eV as shown in Figs. 3(e)—3(h), al¬ 
though both |r^^^(gr)| and |r^"°^^(gr)| take the maximum 
values at qr = 0. These results mean that the susceptibili¬ 
ties should be calculated to obtain the real density-wave wave 
vector in the RG study. 

Appendix C: ANALYTICAL EXPRESSIONS OF THE 
MAKI-THOMPSON AND ASLAMAZOV-LARKIN VERTEX 
CORRECTIONS 

In this section, we show the analytical expressions of 
the leading-order Aslamazov-Larkin VC and Maki- 

Thompson VC for the charge channel, including the 

form factor fa{k). The Maki-Thompson VC is given by 


Xspin ^charge 1 ^= ^d-orb^ represent the (i-orbital spin and 
charge susceptibilities, given by Eqs. (2) and (8), respectively. 
The Aslamazov-Larkin VC is given by 


X^^/iq) = TUjJ2A^{q;p) 

P 


spin 

2^ 




(.p) 


-.X 


charge 




where p = (p, iuj^) and = 2iTnT. Here A/ and AJ are the 
three-point vertices: 


A„(g;p) = Ty]/„(fe+|) 

k 

X {k + q) Gd,d{k - p) Gd,a (k), (C3) 

A'f,{q-,p) = Tj2fp{k' + ^'^Gp,d{k')[GdAk' -P) 

k' 

+ + P + ^)] + q')- (C4) 


Their diagrammatic expressions are shown in Fig. 13. 

For the (i-orbital BDW with (i-symmetry form factor, the 
form factors fa{k) (a = (i, x, y) are set to 


= T^Uj y] /„ (fe + I) + q) GdAk) 

k,k' 

X GpAk') GdAk' + q) fp (V + f) , (Cl) 


where k = (fe,e;), q = {qd^m), H = nT{2l + 1), = 

27rmT, and a, fi = d,x,y represent the orbital indices. Here, 


fd{k) = cos{kP} - cos{ky), fPk) = fAk) = 0. (C5) 

The form factors for the p-ODW with d symmetry are 

fd{k) = 0, /,(/c) = +l, fy{k) = -l, (C6) 

and those with s' symmetry are 

fd{k) = 0, Uk) = ^l, fy{k) = ^l. (Cl) 
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